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A. 

LOCATION-ESTIMATING, NULL STEERING (LENS) 
ALGORITHM FOR ADAPTIVE ARRAY PROCESSING 

CROSS REFERENCE TO RELATED APPLICATIONS 

This application claims priority from the provisional application designated serial number 
5 60/100,291 filed September 14, 1998 and entitled "Location-Estimating, Null Steering (LENS) 
Algorithm for Adaptive Array Processing". 

GOVERNMENT RIGHTS 

This invention was made with government support under Contract NIH-5R01-DC00270 by the 
;e 1 0 National Institute of Health. The government may have certam rights in the invention. 

\| BACKGROUND OF THE INVENTION 

The present invention relates to the field of digital signal processing, and in particular to the field 
|1 to adaptive digital filters for beamfarming transducer arrays. 

Many applications (e.g., radar, sonar, satellite communications, background noise reducing 
hearing aids, etc.) often include a directionally sensitive antenna-array or microphone array 
(generically a beamforming traducer array). These beamforming transducer arrays and their 
associated signal processing are capable of enhancing a desired target source relative to another 
20 non-target source (e.g., a jammer/noise source). The enhancements improve die perceived target- 
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to-janraier power ratio (TJR), which, in turn, improves the performance of subsequent, target- 
based processing. 

One application of a beamforming transducer array is for hearing aids. Hearing aid wearers often 
5 complain about the corrupting influence of background noise upon their ability to understand a 
desired target source. Enhancing the target relative to the background noise reduces this problem. 
Similarly, applications such as satellite commimications, radar and sonar also benefit from the 
directional filtering ability of beamforming transducer arrays. We shall now discuss the basic 
beamforming problem. 

iO 

2 Beamforming Background 

2.1 Basic Beamforming Problem 

1 5 The overall beamforming goal is to extract a desired target signal from within a field of 
non-target, jammer signals using an array of spatially-separated sensors. Specifically, 
beamformers exploit input differences arising from different source directions of arrival 
to create a 'spatial-filter' that preserves the target source while nulling the jammer 
sources. Figure 1 shows the basic beamforming system structure, which filters 
(convolves) each array input x,-[n], f = {0, 1, . . . , M - 1} with a corresponding impulse 
response Wi[n] and sums the resulting signals to form the overall output y[n]: 

M-l 

y[n] = J2 Wi[n]*Xi[n]. (2.1) 

t=0 
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Converting Equation 2.1 into the frequency domain and introducing vector notation 
allows for the beamforming operation to be expressed as an inner product: 

V(u;) = J2 ^'MA'.(w) = K"ii-m^), (2.2) 

1=0 

where 5- A'.(u;), u;.[n] 5 W-{u,'), y[n] ^ Yiu), X(a;) = [A'o(a;), . . . , AV-i(w)]^, 
^ j£(u;) = {Wq{(jj), . . . ,WM-iii^)]^ , and ^ is the hermitian (complex-conjugate trans- 

pose) operation.^ The main issue involved in the design of a beamforming system 
is the selection of the array weights Wlu) that yield the desired spatial-filtering be- 
havior. As described in Section 2.3, these weights may be either time-invariant or 
time-varying. Time-invariant weights typically achieve simple and robust processing, 
while time-varying weights are often capable of superior performance. 

2.2 Beamforming Terminology 

Before discussing the beamforming weight-selection process, this section introduces 
the basic terminology that is used to discuss beamforming systems. 

2.2.1 Coordinate System 

Beamforming systems use spatially-separated sensors to extract a target source from 
within a field of jammer sources. The spatial dimension of this task requires a coor- 
dinate system to describe both array element and source locations. 

Figure 2 shows the three-dimensional, polar-coordinate triple, {po, Bq, 4>o), that 
corresponds to the cartesian coordinate triple (xo, yo, zo). In this representation, (1) po 
is the distance from (ro,yo,^o) to the origin, (2) 9o € [-f, f] is the angle between 
25 i^o, Vo, ~o) and the y^-plane (positive if xq > 0 and negative if xo < 0), and (3) 4>o € 

[--,-] is the angle between the projection (0,3/0,^0) of (xo,2/o,-o) into the yr-plane 
and the positive 2-a.xis. The triples (xo, 2/0,^^0) and ipo.0oM are related follows: 

Po = i/xo + yo + 2^ 

= (2.3) 
<Po = tan-^ ^ 

-^0 
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5 This polar coordinate definition is sligiitly different from traditional polar-coordinate 

definitions, but it simplifies the analysis of broadside arrays (Section 2.2,2), which are 
used for much of this research. 

In the analysis of beamformer behavior, it is assumed that all sources are far-field 
and produce plane waves that propagate across the arra^/, i.e., p = co, which means 
10 that (5, <f>) alone specifies source locations. The theoretical target source location is 

assumed to be known and fixed along the positive r-axis: (5r»^r) = (O^O)- 

2.2.2 Array Configurations 

The array elements lie at known, fixed locations within the array environment and 
may be arranged in any desired configuration. One special type of array is a uniform 
linear array (ULA), in which the array elements are co-linear and equally-spaced by 
a distance d. Figure 3 shows two ULA examples, which are defined relative to a 
specified target location (along the positive z-axis). For these two configurations, 

20 endfire and broadside, the array elements are respectively oriented coaxial with and 
perpendicular to the target propagation direction. 

The simple and regular structure of ULAs allows for straightforward system anal- ^ 
ysis, and, for this reason, the specific systems that are designed and simulated in the 
following chapters will have broadside ULA configurations. The coordinate system of 

25 Section 2.2.1 is convenient for the broadside ULAs of Figure 3b, because only the 
9 location parameter is needed to describe the impact of a specific source upon the 
free-field array (due to the cylindrically-symmetric behavior about the array axis). 
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2.2.3 Source-to- Array Transfer Function, S(a;, 6, 4>) 



10 



The array element locations, the source locations, and the knowledge of the source 
propagation rules within the array environment allow for the definition of a source- 
to-array transfer function that relates each source to each array element. Consider a 
source, 0, 4>), located at (6, (f). This source produces a signal Hi{u, 6, <f>)S{u, 9, 4>) 
at array element i, where Hi{uj,e,4>) is the source-to-element transfer function. The 
overall array input due to S{u>, 9, (f>) is thus 





Xoiuj) 












_ Aa/_i{u;) _ 





Hoiu,9,cp) 
H^^-l(u:,9, <p) 



S{u;,9,4>l 



(2.4) 



where H_{uj. 9, (p) is the source-to-array transfer function vector. Because beamform- 
ing exploits the relative differences between the input signals, the source-to-array 
transfer functions are usually normalized by H^(uj,9,(p). Defining Hi{u:,9,(/>) = 
H,{u},9.o)/ Ho{uj,9,(j>), the normalized source-to-axray transfer function vector is 



Ho{uj,9,cl>) 



1 



(2.5) 



15 Under the assumption of far-field sources that propagate as lossless plane-waves 
across the array and that the array is operating in the free-field_, it is possible to 
determine H_{uj, 9, (}>) explicitly. In this case, the individual-element transfer functions 
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Hi{u!,9, cp) are pure delays, since the amplitude does not decay, 

The delay r,(^,0) is equal to the additional distance, that the propagating 

wave travels between elements 0 and z, divided by the propagation velocity c: 

The array geometry provides a means for determining the di{0, <j}) for a given source 
location. Note that, since element 0 is used as the reference, do{uj, 9, cf)) = 0, which 
yields Ho{uj,9,4>) = 1 as desired in Equation 2.5. 

For a free-field ULA, each di{e,4>) may be stated in terms oidi{6,<l>). Specifically, 
the equal inter-element spacing and the linear array orientation result in 



di{uj, e, 4) - di.iiij, 9, <f>) = di{uj, 9,<f>), i = {1, 2, . . . , M - 1}. 
This reveals that di{e,<p) = idi{9,(f)), and, since Ti(9,<f>) = ^i^M^ this yields 





1 




1 




1 














Hi{<^,9,4>) 
















■ (2.6) 




.. .(Af-i)<f,(e.*) 













Based on the coordinate system in Figure 2 and the ULA examples in Figure 3 
the following di{6,<f)) values can be determined, 



Endfire: dUO, <f>) = d cos6 cosS, 

(2.7) 

Broadside: di{9, <i>) ~ d sinO, 



where d is the inter-element spacing. 
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2.2.4 Overall Input Signal 



The source-to-array transfer functions allow for the expression of the overall array 
input as the superposition of the individual-source array inputs. For example, given 
L sources 5(u?, 01, c^/), / = {1, 2, . . . , L}, and the corresponding source-to-array transfer 
functions H{uj,9i^(pi)^ the array input is 

Since the following beamforming techniques are narrowband and operate only on a 
narrow band of frequencies, the following discussions suppr^^ss the explicit u depen- 
dence. 



2.3 Beamforming Methods 

In order to see how W determines a beamformer's spatial filter characteristic, combine 
Equations 2.2 and 2.8 



^t.S{Ou<i>i)W^H{eu<f>i)' (2.9) 



This output expression reveals that each source S{9u(f>i) in the array input is also 
present in the array output, but it is scaled by the location-dependent factor known as 
the directional array response G{d, (j>) = W^K{&^ 4)^ This scaling factor demonstrates 
the effect of E£ upon a source from the 4>) direction. 

The overall spatial-filtering, beamforming goal is to determine a ]£ such that 
G[9, 4>) = I for the target source and G{9, cj>) = Q for the non-target, jammer sources. 
In the special case where the target and jammer source-to-array transfer functions 
are known [i.e., the target is at (^Ti^t) ^-^d the jammers are at known locations 
{9j,k^4>j,k)'} ^ = {1,2, ...,M — 1}], it is possible to formulate a system of linear 



equations whose solution achieves this goal: 
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G{9T,<i>T) 



H^{9t, <pt) 




1 




W = 


0 






0 



(2.10) 



where the equality G{e, 4>) = W" H{d, <f>) = H{9, <i>)"W has been used. 

Equation 2.10 cannot be used for general beamforming purposes, since the jam- 
mer locations are unknown and, usually, time-varying. Other methods have been 
de\'eloped to generate beamforming weight vectors that meet the desired goal. These 
techniques fall into two categories: 

^- Time-Invariant Systems: in which W_ is fixed and optimizes a general perfor- 
mance criterion, such as the array directivity index. 

Time- Varying Systems: in which ]£ is continually updated to suit the particular 
beamforming environment. 

Time-invariant systems have the ad\'antages of simple implementation and robust op- 
eration, but the unchanging nature of these systems requires that the beamforming 
be optimized to yield generally-good performance in all potential array environments. 
This requirement prevents these beamformers from achieving the maximum possi- 
ble level of jammer cancellation that is attainable for the particular environment in 
which the system is operating. Time-varying systems, on the other hand, have the 
advantage of continuously updating the beamforming in order to achieve the optimal 
performance for the particular array environment. This improved performance arises 
at the cost of more computationally-demanding processing, and these systems may 
cancel the target under certain conditions. 
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Since LENS processing is time-varying, this research considers only time-varying 
systems. The following discussion introduces these systems in three parts. First it 
presents the simplest concept of optimal beamforming for a given environment and 
introduces the robustness problems that arise with this system. It continues with 
5 a discussion of the traditional, robustness-enhancing modifications to the original, 

non-robust weight-vector definition. It concludes with a description of some imple- 
mentations that may be used to approximate these systems. 



2.3.1 Traditional, Non-Robust Beamforming 

10 Equation 2.11 describes the traditional, non-robust frequency-domain beamforming 

weight-vector^ which is chosen to minimize the expected output power at each fre- 
quency, subject to the constraint of preserving unit gain for the known target- 
source location: 

IZXR = argmin E[\Y\^] = argmin E[W^XX^W] = argmin W^RxxK 

subject to (2.11) 

15 where Rxx = E[ XX "] Is the input covariance matrix and target — Ki&T, <t>T) Is the 

target-to-array transfer function vector. Solving this problem via Lagrange multipliers 
yields the solution: 

ii^arget A' X iiiarget 

In beamforming literature, this weight vector is referred to by several different names, 
20 including the Maximum Likelihood Method (MLM) weight vector and the Minumum- 

Variance, Distortionless Response (MVDR) weight vector. This research uses a third 
appeiation, NR, which reflects that this is the baseline, non-robust situation of the 
Constrained, Minimum-Variance (CMV) weight vector [described in Section 2.3.2]. 



10 



This definition of the optimal beamforming weight vector seems logical - by mini- 
mizing the expected output power while preserving the target signal, the weight vector 
must be attenuating the jammer sources in the environment. This system works well 
under ideal conditions, but, unfortunately, problems arise under more realistic con- 

5 ditions. These problems result from discrepancies between the assumed and actual 

target-to-array transfer function vectors: target target- 

The primary effect of the discrepancies between fl^arget ^'^^ target is that they 
can result in cancellation of the target signal - an obviously undesirable effect in any 
beamforming system. To see how this behavior arises, consider the case when only 

10 the target is active within the array environment. If H^arget -target different, 

then it is possible for the beamformer to find ]£nr such that 

^S^^rgetli^NR = U 

:l ^rg.t&R = 0- 

% In other words, the beamformer can completely eliminate the actual target signal 

M while still preserving the assumed target-to-array transfer function. 

Generalizing this result to situations in which both target and jammer sources 
are present, the following behavior is observed. When the target is weak relative to 
S the jammers, i.e., the TJR (target-to-jammer power ratio) is low, then the jammer 

I power dominates the beamforming, and the system cancels the jammer signals. When 

^ the target is strong relative to the jammers, i.e., the TJR is high, then the target 

20 dominates the beamforming and the target cancellation described above may happen. 

The most obvious means of avoiding target cancellation would be to obtain a better 
target-to-array transfer function estimate [i.e., to reduce the discrepancies between 
target ^^rget] f^r usc iu Equatiou 2.11. Unfortunately, this is not an easy task, 
given the following causes of discrepancies between target target- 

25 1- Array Mis- Alignment/Mis-Steer: in which the target is not located exactly at 

its assumed location, e.g., it might be near-field or slightly off center. 

2. Element Placement Errors: in which the arrav elements are not located at ex- 
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actly the specified array locations. 



^- Element Gain Misma.tr^r in which the array elements have unequal gains. 

Non-Free-F ield Array Operation: in which objects that are located near the 
array scatter the incoming wave-fronts, e.g., a hearing-aid array mounted on 
the head. 

^- Multipath-Propagation/Reverberation: in which the target signal arrives from 
both the actual target location and non-target locations (wall reflections). 

These effects are difficult to estimate and usually change over time, which renders 
it impossible to maintain a highly accurate estimate of EU^,^,, in practice. For this 
reason, modifications have been made to the definition of the optimal beamforming 
weight vector of Equation 2.11 in order to increase its overall robustness to target- 
cancellation. The following section describes one class of modifications. 

2.3.2 Robust Processing: CMV Beamforming 

Constrained, Minimum- Variance (CMV) beamforming is a classic approach to in- 
creasing beamforming robustness. It operates by imposing additional weight-vector 
constraints to the traditional, non-robust optimization of Equation 2.11: 

Kcm' = argmin W^RxxEL ('> 13) 

subject to: 

C"W = g and W^W < S'^, 

where (C,fir) describes a set of linear constraints (that must include preservation 
of the assumed target-to-array transfer function, HS,,^,,W = 1) upon the weight 
vector and S'^ is a quadratic constraint on the weight vector norm. Note that CMV 
beamforming reduces to traditional, non-robust beamforming [Equation 2.3.1] when 
C = .target, £ = 1, and = oo. Additionally, when 5-^ = oo, then no quadratic 
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constraint is Imposed upon the weight vector and CMV beaniforming reduces to 
linearly-constrained, minimum- variance (LCMV) beamforming. 
The solution to Equation 2.13 may be vvritten as: 



Kcmr = {Rxx + dr'C[C''{Rxx + dr'C]-'g. (2.14) 

When there is no quadratic constraint, i.e., when S'^ = oo such as for traditional, 
non-robust and LCMV beamforming as described above, this solution is closed-form 
with e = 0. When the quadratic constraint is active, however, this solution is not 
closed-form and e must be adjusted iteratively until W^W < 5'^. 

The constraints C^W = g and W^W < S"'^ discourage target cancellation in 
the following manner. The linear constraint C^Wi — g imposes restrictions on the 
behavior of the directional array gain G{6, <f>). For example, the required constraint of 
afar^et^ — 1 ensures that the assumed target gain is unity: G^(^ri ^t) = 1- Derivative 
constraints may be added to force G{0^<(>) to have a flat derivative about the 
assumed target location, which prevents the S5''stem from canceling targets that are 
close to the assumed target location. Alternatively, eigenvector constraints that keep 
G[6,<f) 'near' a desired behavior about the assumed target location may be used, ' In 
general, the linear constraints are designed to increase system robustness to highly 
structured deviations (e.g., array misalignment) between the actual and assumed target- 
to-array transfer function. 

The norm constraint, IV'^iy < 5"^^ imposes restrictions on the white-noise-gain 
of the system. Specifically, the white noise gain, measures the 

improveincnt in the output relative to the input TJR (target-to-jammer power ratio) 
when no <iirectionai jammers are present and only uncorreiated noise is introduced 
at each array element: 

_\eL^_ 1 

where the fact that ^^^getB^ = 1 has been used. By maintaining W_"W_ < 5~'^, 
then Gy, > 6^, which guarantees a minimum level of beamforming improvement, even 
in the presence of uncorreiated sensor noise. This norm constraint increases svstem 
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robustness to factors such as array-element placement error, which are independent 
and may be regarded as additive uncorrelated noise at the array input. 



2*3.3 Traditional Implementations 

In general, the array environment and input covariance matrix are unknown and 
time-varying, which prevents direct solution for the optimal CMV weight vector as 
described in Equation 2.13. Approaches have been designed, however, that approx- 
imate the optimal solution and that track the changing array environment. This 
section introduces three such implementations. The first two implementations are 
possible only when = oo, such as for traditional non-robust and LCMV beam- 
forming sub-classes of CMV beamforming described above. The third technique is 
more general and may be used with any CMV formulation. 

All implementations use a temporal stream of array-input observations, to 
derive a corresponding temporal stream of beamforming weight vectors, {W^^^y}. 
For the broadband processing simulated in much of this research, the narrowband 
^^^^ arise through block FFT processing on the broadband input. 



Direct Matrix Inversion: In the absence of a weight vector norm constraint, i.e., 
when (S""^ = co, the solution of Equation 2.14 is closed form. In this case, EL^y^iy 
may be obtained by forming an input-covariance estimate Rj/j^ 

and using this estimate in Equation 2.14. This method is known as the 
Direct Alatrix Inversion or DMI method, because it involves the inversion of the esti- 
mated covariance matrix. The main advantage of this implementation is its relatively 
fast adaptation time to a changing array environment: adaptation time is governed 
entirely by the delays associated with the covariance matrix estimation process and 
are independent of the array environment. Unfortunately, this fast adaptation comes 
at the cost of matrix inversions, which can be computationally expensive [0{M^) 
operations] . 
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Frost Processing: An alternative impiementation to the DMI 

implementation, also for the special case when 5'^ = 00, is a Least-Mean- 

Squared (LMS) based adaptive approach that uses a gradient search to converge upon 
the optimal beamforming weight vector. Specifically, the weight vector follows the 
5 gradient towards minimal output power, but it is projected back into the constraint 

subspace (defined by C^W = g) at each iteration: 

where 

\JU = C{C^Cy^g and ^ I ~ C{C^C)''C^ . 

This processing has the advantage of being less complex [0{M^) operations] to imple- 
ment than DML but this complexity advantage comes at the cost of longer adaptation 
10 times: adaptation time is proportional to where Xmin is the smallest eigenvalue of 

the underlying input Rxx^ which means that certain array environments can exhibit 
very long adaptation times. Details of this technique are set forth in the paper 

entitled "An Algorithm for Linearly Constrained Adaptive Array Processing", Proc. 

IEEE, 60:926-935, August, 1972, O.L. Frost. 

15 Scaled Projection Algorithm: The processing set forth above may be altered yield 
an adaptive beamfcrmer implementation that includes the weight vector norm con- 
straint, W^W < <f~^. In this case, the basic Frost algorithm is used, but the weight 
vector is restricted to obey W^W < 5^'^ at each step: 

Step 1: £^'^^' = [VI^O _ tiX}'^2L^'^»w!^^^] 

Step 3: = + Zi+i (2-1') 



where: = C{C" C)-'g, Pc = I - C{C^C)-'C, 
and b'' = 5-^ -g"(C^C)-'g. 
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As with Frost processing, this processing has the advantage of being relatively simpi 
to implement [0{M^) operations] but it may exhibit long adaptation times. 

As set forth above, a significant about of care (i.e., constraints) must be taken when 
defining the optimal beamforming weight vector in order to yield a solution that is 
robust to situations that may cause target cancellation. The CMV definition provides 
one possible solution to this problem, but it suffers from some disadvantages. 

1. The CMV constraints are indirect - their unpact upon the beamforming 
performance and their prevention of target-cancellation can be difficult to 
understand. This renders the overall constraint-selection process 
confusing. 

2, Each luiear CMV constraint in {C,g) reduces the beamformer degrees 

of freedom by one. Under good beamformmg conditions, e.g., when the 
TJR is low, this reduction in the degrees of freedom limits the potential 
jammer cancellation. 

Therefore, there is a need for a robustness-control mechanism, which yields a system 
that has direct and flexible control over the beamforming process. 
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SUMMARY OF THE INVENTION 

Briefly, according to the present invention, an adaptive multi-tap frequency domain digital 
filter processes an input signal vector X from a plurality of spatially separated transducers that 
detect energy from a plurality of sources including a target energy source and at least one non- 
target energy source. The filter receives and processes the input signal vector X to attenuate noise 
from non-target sources and provides an output signal vector Y. Tap weights for filter are 
selected by first parameterizing each of the tap weights such that each of the tap weights 
is characterized by a vector of parameters P^^^ , and the solving for each parameter of the vector 
^opt^y minimizing the expected power of the array output signal Y. A robusmess-control 
transformation is then applied to the vector B to provide a robust vector B , wherein the 

—opt Urob 

robustaess-control transformation identifies and reduces target canceling components of the vector 
Lpt"^^^^ preserving non-target canceling components. Finally, the weight vector indicative of 

the filter tap weights is formed as a function of the vector B . 

i—rob 

The tap weight selection technique of the present invention preferably executes in real- 
time. The technique may be used with ttansducer arrays such as microphones (e.g., hearing aides, 
sonar processing, etc..) and radar receivers where it is desirable to attenuate signal contributions 
associated with non-target signal sources. 

Beam forming parameterization identifies parameters of the vector B that describe the 

—op/ 

nuUmg behavior of the tap weights . 

The present invention uses a direct and flexible robusmess control mechanism to yield a 
beamformer that provides good performance and it robust to a wide variety of adverse conditions. 
Notably, the present invention includes a novel robustoess-conirol mechanism, which 
parameterizes the beamforming problem and yields a system that has direct and more flexible 
control over the beamforming process. 

Advantageously, selectmg the tap weight of the adaptive filter according to the present 
mvention nulls (i.e., attenuates) signal components associated with non-target noise sources, while 
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retaining signal components associated with a desired target source(s). 

These and other objects, features and advantages of the present invention will become more 
apparent in light of the following detailed description of preferred embodiments thereof, as 
illustrated in the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram illustration of a beamforming array; 
Figure 2 is a plot of a coordinate system; 

Figure 3A and 3B (collectively Figure 3) are simplified pictorial illustrations of two 
uniform linear array (ULA) configurations with respect to a target; 

Figure 4 is a block diagram interpretation of the processing according to tiie present 
invention for a four element ULA; 

Figure 5 is a plot of the transformation of p^^, , to ^ ; 

Figure 6 is a pictorial illustration of the distance metric Q between the vector Z and the 
subspace S; 

Figure 7 is a plot of the contours of various values of Q,{0) for ULA arrays having 2, 3, 
4 and 5 elements; 

Figure 8 is a plot that illustrates three LENS parameter classification regions for an array 
having three elements as described by Equation 3.20 with =30° and = 60°; 

Figure 9 is a plot of the contours for a three element array and Q{fi^ ) =30°, 40°, 50° and 
60° with indicated as (x); 

Figure 10 is a plot that illustrates three LENS parameter classification regions for a three 
element array according to Equation 3.23 when and are chosen to reflect the R^. =30° and 
Rj =60° contours, respectively; 

Figure 11 is a plot of illustrating four examples of fi^^,^ for a two-element array along with 
the corresponding LENS robustness regions; 

Figure 12 is a flow chart illustration of processmg steps accordmg to the present invention; 

Figure 13 is a flow chart illustration of alternative embodunent processmg steps accordmg 
to the present invention; 

Figure 14 is a flow chart illustration of yet another series of alternative embodiment 
processing steps according to the present invention; 
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Figure 15 is a block diagram illustration of processing steps for solving for yS^^ given a 

known covariance matrix; 

Figures 16A and 16b (collectively Figure 16) are block diagram illustrations of processing 
steps regarding the adaptive relaxation based solutions; and 
5 Figure 17 is a plot that illustrates the relative computational complexity for the complete 

update relaxation, partial update relaxation and gradient LENS implementations. 
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DETAILED DESCRIPTION OF THE INVENTION 



3 LENS Beamforming Processing 

The present invention is directed to a beamforming technique that modifies traditional non-robust 
beamforming to increase beamforming robustness by preventing target cancellation. Significantly, 
unlike the prior art, the present invention solves for the robust weight vector indirectly by dividing 
the beamforming process mto three components/steps. 

Like CMV beamforming [Section 2.3.2], Location-Estimating, Null-Steering (LENS) 
beamforming modifies traditional, non-robust beamforming [Section 2.3.1] in an at- 
tempt to increase beamforming robustness by preventing target cancellation. LENS 
differs from CMV processing however, in that it does not solve for a robust beamform- 
ing weight vector directly. Instead, it solves for the robust weight vector indirectly 
by dividing the beamforming process into three components: 

L Estimation of the optimal LENS parameter set, that describes the tradi- 
tional, non-robust beamforming problem. 

2. Transformation of into a 'robust' parameter set, that is designed to 
ensure robust beamforming and to prevent target cancellation. 

3. Use of ^^^^ to generate a final beamforming weight vector. 

It should be noted that LENS processing assumes that target = [1? 1? • • • ^ 1]^ = i- 
This target-to-array transfer functions arises naturally for broadside ULAs [Sec- 
tion 2.2.2]. For more general arrays, it may be realized by scaling each array input i 
by the steering factor jg^^^^i where {target}** is the i^^ element of fftarget- 

The primary advantage of the LENS parameterization is that it allows for the use 
of an innovative and direct method of beamforming control. Consider as an example 
Uniform Linear Arrays (ULAs), for which each parameter /3opt,i in determines a 
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beamforming null location in the array environment [Section 3.3-1]. In the process 
of transforming into it is possible to control where the final beamforming 
weight vector will place its nulls [Section 3,3.2] and to prevent target cancellation 
by restricting these nulls to be a certain angular distance away from the assumed 
target location. (It is this manifestation of the algorithm that leads to the Location- 
Estimating, Null-Steering appellation.) 

This chapter develops the basic components of LENS beamforming. It begins 
by introducing the the LENS weight vector and parameter set. It then determines 
the optimum LENS parameter set to describe traditional, non-robust beamforming 
and defines the LENS robustness transformation. It concludes by drawing these 
components together into the overall, general LENS algorithm. 

3.1 The LENS Beamforming 
Weight Vector 

The key component of LENS processing is the LENS beamfonning weight vector. 
IilLENS,A/(^)- il/-element array, this weight vector is determined from the 

M - 1 LENS parameters, P = [/3i,/?2, • • . ,/?a/-i], as follows: 

Ii:LENS.Af(^)'' = {^J' \{l-j/3i) (l+i/?i)]* (3.1) 
[(J-j/?2) (l+i/?2) ]*•••*[ (l-j/?M-l) (l+i/?A/-l) 

where the '*' operation denotes row vector convolution. It turns out that the LENS 
parameters ^ provide a convenient description of the beamforming that is realized by 
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BLlE^Sm(§) ^^^'> consequently, they may be used to control the beaniforming and 
to prevent target cancellation. 



10 



15 



In order to understand how W^lENSM^S) ^^^^^^ obtain a 'first look' at the 

relationship between 0 and EIlENS^uC/')? consider the two-element W LENS,?(^0' 



(3.2) 



Defining ib such that 



0 = 



01 'j 



(3.3) 



where both 0 and Pi may be complex, and substituting into Ii[LENS,2(/^i) Y^^Ws 



I£lens,2(.^0'' = 



1-0 



where iELENS.2(A)'%ar6et = 1 and I£lenS.2(/?i)'^^ = 0, with 



1 

The one-to-one relationship between 3i and shows that /?i completely describes the 
nulling behavior of W.LENSm^§)' 
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The precise nulling behavior depends upon the actual source-to-array transfer 
functions HS9^<j>) for the array in question. Consider the special case of far-field 
sources, for which the source-to-array transfer functions are given by Equation 2.6: 



1 



By equating il? = e-''^'"^^'*^) and substituting into Equation 3-3, then ELLENS, iV/(^) 
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clearly nulls source locations [9, (j>) that satisfy: 

In this special case, since T{d^<f>) is real, /3i is also real. As discussed in Section 3.3.1, 
for general .IZ-element ULA architectures with far-field sources, each LENS parame- 
ter /?,■ determines one beamforming null location through a similar relation to Equa- 
tion 3.4. For more general non-ULA array architectures and for near-field sources, 
the relationship is more complex. 

The iterated-convolution, M-element LENS weight vector, ]ilLENS,iV/(/^)? 
interpreted as a cascade of the two-element LENS weight vectors [Equation 3.2]. 
Figure 4 shows this cascade for the special case of a ULA array architecture with 
far-field sources, for which each LENS parameter and; consequently, each stage in the 
cascade produces one null in the array directional response [Section 3.3.1]. 



3, 1.1 LENS Weight Vector: a Re-Formulation 

The LENS beamforming weight vector 1ElENS,m(^) re-expressed in a man- 

ner that simplifies the analysis of the LENS algorithm. It can be shown that 

IilLENS,A/(^) equal to the linear transformation Tm of another vector £iLENS,A/(^)- 

1£lens,m(^)'' = UjLmsMiifT^^ (3.5) 

where £Z!lENS,m(^) defined as 

ULEmM{§)'' = (2 j t 1 A ] * [ 1 /?2 ]*■••*[ 1 /?A/-i I (3-6) 
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and T\{ is given by 



[11] * [11] * 

[11] * [11] * 

[11] * [ - j j ] * 

[ -i J ] * [ -j j ] * 



* [11] * [11] 

* [ 1 1 ] * [-3 j\ 

* i-j j] * I -j j ] 

* i-j j] * [ -j j ] 



k — I pairs per row 



3.2 Optimal Beamforming Re-Expressed as LENS 

The LENS weight vector allows for the LENS parameterization of traditional, non- 
robust beamforming [Section 2.3.1]. Specifically, define to minimize the expected 
power in the LENS array output signal Y = W lenS,jU (^)^^' 

3 = argmin E[\Ym'] = argmin Kj^ENSMii)'' RxxKlENsMi)- 

This yields a ii£LENS,A/(^^pf) ^^^^ equal to the traditional, non-robust IHnr, 
from Equation 2.11. To see why this is so, consider the following proof. 



Given that satisfies Equation 3.8, then WjLENSM(§^pt) ~ ^NR- where 
li^NR satisfies Equation 2.11 in the case where steering delays have been applied 
to the array so that ^target = 1- 

KnR = a^rgmin W^RxxK (3.9) 
subject to: 1^1£= L 
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Define the set A of M-dimensional beamforming weight vectors that are rep- 
resented by ]£LENS,Ar(^) as ^ spans the space of M - 1-dimensional complex 
vectors, C'^^-^: 

-4 = {IZlens.a/(/?) € C^'l^ 6 C'^^-^}. 

It can be shown that A is equivalent to the set B of ail iW-dimensional 
complex vectors W that sum to 1: 

where i = [1, 1, . . . , 1]^. Hence A = B. 

Using the set B, Equation 3.9 may be re-written as: 

Knr = argmin W^RxxK, 

where W clearly varies over the entire set B for the optimization. Repeating 
Equation 3.S, the optimal LENS parameter set is: 

(Lpt = aJ^gra^in IElens.a/(/?)'^/2.va-I£lens.a/(.^)' 

where ]£leNS..\/(,^) varies over the entire set A for the optimization. 

Since (1) W varies aver B, (2) KlenSmH) varies over A, and (3) A = B, then 

Ii:LENS,M(^„,) = W^R. 



3.2.1 Closed-Form LENS Solution 

The equality of ]£lens..v/(^) and W-^^, as described above, and the LENS re- 
formulation from Section 3.1.1 may be combined to yield a closed-form expression 
icpf Specifically, using the closed-form solution for 1£^TR [Equation 2.12] and 
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recalling that fl^arget = L then 

Now, recall that i:LENS.A/(£,pJ^ = ULENS,Mitpt)"'^M [Equation 3.5] and 
it is known that Qm = Tj^j- exists, which leads to 



X X ~ 



Finally, use the fact that the elements of the vector £^lENS,m(^)^ iovm the coefficients 
of a polynomial that has roots at — [Appendix B] to obtain the expression for 

= -roots {U^EmMii^f) = -roots (^^^) • (3-10) 

Equation 3.10 reveals that may be obtained by inverting Rxx and performing 
a polynomial root solution on the polynomial described by I^lENSmC/')- 



3.3 Transforming into ^^^^ 

As shown above, is a parameterization of the traditional, non-robust beam- 
forming of Section 2.3.1. This means that a beamformer that uses WjiENSM^§^pi^ 
would exhibit all of the non-robust, target-canceling characteristics of the original, 
non-robust approach. The key feature that differentiates LENS beamforming from 
traditional, non-robust beamforming is the robustness-control transformation of 
into This transformation uses the relationship between ^ and the nulling be- 

havior of W£lENS,m(^) to identify and eliminate any target-canceling components of 
§^pi ^^^^^^ preserving the jammer-canceling components. 

The robustness transformation is achieved by classifying each optimal LENS pa- 
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rameter (3opt^i into one of three sets and then transforming it according to the rule: 



fSopt, € J (3.11) 



where T and J are the 'target-canceling' and 'jammer-canceling' LENS-parameter 
regions, respectively. Figure 5 demonstrates this process graphically for param- 
eters in each of the three classes. For /^opta 6 T, the LENS parameter contains 
primarily target-canceling information, and so /?ro6,i is set to 0, which is a 'safe' 
non-target-canceling parameter value. For popta € the LENS parameter contains 
primarily jammer-canceling information, and so /?ro6,i is set to the unaltered 



For Popt,3 £ T li J ^ the LENS parameter is in an intermediate region and (3rob,i is 
10 set to F[/3opt,i], which is a smooth transition between the ^opt.i € T and (3opt,i € J 
extremes. 
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Equation 3.11 reveals that the main element in the transformation is 

the parameter classification T, JT, and TU J. The following discussion develops this 
classification in two stages. The first stage explores LENS parameters and investi- 
gates the beamforming information contained in while the second stage uses this 
information to derive the actual parameter classes. These two components are then 
combined into the final LENS robustness transformation, and some simple example 
situations are explored to demonstrate its operation. 

3.3-1 LENS Parameter Classification, Part 1: 
Beamforming significance of ^ 

In order to classify a LENS parameter as either target-canceling or jammer-canceling, 
it is necessary to understand the significance of /? as. it relates to the overall system 
beamforming behavior. This relationship is stated as a fact for the two-element array 
of Section 3.1 [Equation 3.3], but it is more complex for general, M-element array 
architectures. 

For general axray architectures, P and ]£^leNS,m(/^) related by the fact that 
^ describes a basis for the null-space, A/'[I£lENS,a/(^)]9 WjleNSm(§)-> ^^^^^^ 
turn describes the source-to-array transfer functions that are nulled by 1ElENS,a/(^)- 
Given ^ = [^i, /?2i * • . , 0m^i]'^^ this basis consists of the vectors 



, i = {l,2,,..,M-"l}. (3.12) 



To see why {Y^(§)} forms a basis for J^[W lenS,m(^M- consider the following: 

1. Af[W i^]<\sm(§)] an {M — l)-dimensional space: This arises from the fact that 
the null-space of a single M-element vector is always (M — l)-dimensional. 
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2- span[jA(jg)] is an {M - l)-dimensionaI space: Use the A/ - 1 vectors (^) and 
the vector /^target to form the columns of an M-by-M Vandermonde matrix. As- 
suming that the are unique , then this Vandermonde matrix is non-singular, 
and its columns, including the Vi{0), are linearly-independent . These M - 1, 

5 linearly-independent Vi{0 then span an M - 1-dimensional space, 

3- span[V^(^)] C A/'[IELENS.A/(/g)]- It can be shown that each Vj{§) is 
an element of Ar[lFLENS,A/(^)], i.e., that Ki^^sMiB^^d) = 0- For this 
reason, any vector that may be expressed as the weighted sum of Y^{0) - that 
is to say any vector within span[V£,-(^)] - is also within -Ar[]£LENS,iJ/(^)]- 

10 The only way that one (Af - l)-dimensional space may contain another (M - 1)- 
dimensional space is if the two spaces are equal, and so these three properties reveal 
that span[]/ (^)] = J^IKiy^^smW- This is equivalent to stating that the vectors 
y\i(0) provide a basis for J^[KLEm,M(i)]- 

This leads to the conclusion that the beamforming weight vector W£lENSm(/^) 
15 will null any source-to-array transfer function H16, 4>) for which there exists a set of 
complex coefficients, {c,}, such that: 

M-l 

H{e,<f>)= (3.13) 
1=1 

Since the target-to-array transfer function Sj^arget = i is known, it is possible to use 
the V_,{§) to determine how 'close' target is to Ar[l£LENS,A/(^)] = span[V.(^)] and 
20 to define the desired LENS-classification regions T, J", and TTTj. Before proceeding 
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with this process, however, it is interesting to consider the special case of ULA beam- 
former architectures with far-field sources, which display a convenient correspondence 
between the LENS parameters and the beamformer null locations. 
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5L5 



Special Case: ULA Architectures 

Recall the ULA source-to-array transfer functions for far-field sources 





1 




1 



























(3.14) 



where T{9.(f>) is the additional temporal delay experienced by a source from 
when traveling from array element i to array element {i + 1). It turns out that, for 
this special type of source-to-array transfer function, 1ZlENS,m(^) ^^^^^ only source 
locations with 11(0^ (/>) exactly equal to one of the basis vectors 



(3.15) 



for some i € {1, 2, . . . , M - 1}. 



The <= property of this relation is straightforward to see: since Y^iff) is a basis vector 

for AflKiEmMii)] tlien^j:^(^,<^) =£.(^) is clearly within y[]£LENS.M(^)]- The 
=^ property is less obvious, however. Consider the M-by-M A matrix formed using 
the M -I vectors and the vector HjjiJfi,4>) € ^I]£LENSjvf(£)]- 



njjLA 



{OA) 
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Since {Vi{§)} forms a basis for -VIIElENS.m (^)], then ^i,>i(^, <t>) € ^[WlENS-mC^)! 
must be linearly dependent upon the 1^(£), which means that A is singular. Given 
the structure of Vi{§) [Equation 3.12] and HjjLA{d,<i>) [Equation 3.14]), then /I is a 
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Vandermonde matrix, which is singular only if two of its columns are equal. Since 
the 0i and y^i0) are unique, then the Vandermonde A is singular only if H^uXG-, <p) = 
Y^{§) for some which proves 



By equating the elements of ]£,■{/?) and H_{6, (f>), Equation 3.15 yields 



Hiie, <j>) = e-''"^^^'*) = for some i = 1, 2, . . . , M - 1, 



which may be re-written as: 



j3i = cot — , for some i = 1, 2, . . . , M - 1. 



(3.16) 



(3.17) 



Solving this relationship for each f3i reveals that WlENS,Mi§) ^^^^^ sources with 
locations that satisfv: 



{Oi, ^i) =: r-^ cot-^ a) , i = 1, 2, . , . , M - 1. 



(3.1S) 



Note that these results match those stated in Equation 3.4 for two-element arrays 



Given that minimizes the expected beamformer output [Equation 3.8], then 
I£lENS,.u(j^p^) should null the non-target sources in the array environment. This 
knowledge together with Equations 3.17 and 3.18 suggest that the are, in some 
sense, a set of non- target location estimates. 

It should be noted that the special ULA relationship between the physical null location 
for WiiENSM{fi) LENS parameter set j3 holds for any array architecture that has 

a source-to-array transfer function of the form: 

1 
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Generally, only free-field ULAs will exhibit this form of source-to-array transfer 
function for all source locations. For other arrays, however, it is possible to pre- 
process the array input to yield source-to-array transfer function structures that are 
similar to the ULA structure for a subset of source locations that are of interest for the 
5 particular beamforming application. 



10 



15 



20 



In this case, H{e,cl>) € J^iKLEmMi)] 

H{0,<i>) = Vi{0), i = {l,2,...,M-l}, 
which leads to nulls at locations that satisfy 

^,(,,^) = |±i =^ (,,^^) = ff-g±i), i = {i,2,...,M-i}. 

Note that, by setting Hi{d,<f>) = this result yields Equation 3.3 as stated for the 
two-element array in Section 3.1. 

3.3.2 LENS Parameter Classification, Part 2: 
Generating Parameter Classes 

Given the preceding relationship between §_ and W:LENS.Af(^. i* is now necessary 
to develop the classification procedure that labels LENS parameters as either T, 
J, or T\JJ. This classification is achieved by measuring the impact of /?f upon the 
'closeness' of the assumed ^t^rget to UlW^-EmA&i'' (1) i^^« ^11£lENS,m(^)] 
to be close to ^[iarget, then /?i is considered to be target-canceling, i.e., A 6 T, (2) if 
(3i causes J^{W^-mS,M^§)\ to far from H^^,, then A is considered to be jammer- 
canceling, i.e., A € J, and (3) if A causes Ar[W:LENS.M(£)] to be an intermediate 
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distance ''fom ^arget^ then has an intermediate efFect, i.e., /?,• € T U J . 



Distance Metric, ^\{^) 

One convenient distance measure between ^V[ii:LENS,A/(^)] and target is the angle 
between the subspace and the vector. Figure 6 shows the concept of Q for 
a hypothetical subspace S and a vector Z, where Z has been separated into two 
orthogonal components, Z5 € 5 and Z- Z5 J. <?, that sum to form Z. When Q = 0, 
then Z5 = Z, which means that S contains Z entirely, and the distance is 0. On the 
other hand, when H = 90°, then Z5 = 0, which means that S contains no information 
about Z, and the distance is maximal. 



In practice, the dependence of Q(^) upon §_ is complicated, and it is difficult to 
determine the impact of an individual /?,- upon this metric. It is possible, however, to 
define a set of metrics, where is equal to the angle between target 

! and the span of the single null-space basis vector 

i = {l,2,...,M-l}, 

where | • | is the standard complex scalar magnitude, || • Ija is the Euclidean complex 
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vector norm , and the second equality comes from the knowledge that ^fftarget ~ 1 
and that | [^target 1 12 = V^- Although the set {fti{p)} does not provide as accurate 
a distance measure between ^\f[WiLENSM(§)] target that provided by 
it does provide a general idea of this distance. Additionally, it has the advantage 
5 of being relatively simple and of measuring the effects of each LENS parameter f3i 

individually. Note that, like ^iiP) ranges from a minimum-distance value of 0"" 

to a maximum-distance value of 90^. 

Ideal Parameter Classifications 

10 The metric serves as the vehicle through which 0opt,i is classified. Figure 7 

shows contour plots of fij(^), as /?j varies over the space of complex numbers, for cases 
of M = 2, 3, 4, and 5 array elements. Specific contours indicating = 30^, 60^, 90° 

are shown. Note the following properties of the metric Cti{0 and the information that 
it contains about f^ii 
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1. achieves its maximum value of 90° at precisely Af - 1 points, known as 
the 'orthogonar values of /?,■: p^, ^- = {1, 2, . . . , M - 1}. At these values, the 
span of the individual y_-{^ is orthogonal to targets and /?,- is considered as 
containing jammer-canceling information only. 

5 

2. n.(^) -J- 0° as oo. At these values, span{V;,(^)] practically contains 
^target) ^'^^ Pi is considered as containing target-canceling information only. 

3. riild) decreases regularly as the distance between /?.• and the orthogonal values 
10 increases. At these intermediate values, ^ contains both target-canceling and 

jammer-canceling information, with parameters near to 0^ being more jammer- 
canceling and parameters with large |/?,| being more target-canceling. 

These observations suggest the following classification of 

niiP) < fir =^ A € r, 

15 ^j<ni{§) =^ /?, €jr, (3.20) 

nj<ni{§)<nT /SieTuj. 

Figure 8 indicates the three classification regions for an M = 3 element array when 
Qf = 30° and Qj = 60° , Note that the specific values of Qr and Q, j used in a specific 
LENS implementation will depend upon the array architecture and the frequency of 
the narrowband beamformer input. Boundary selection is discussed in more detail 
20 below for the case of ULAs. 

A note on the 'orthogonal' values of/?,-: It is possible to calculate explicitly the 
M - 1 orthogonal values of /3i that result in = 90°. Specifically, these are the 

/3i that yield = 0. Defining a,- = and recalling the structure of V;-(^) 

25 [Equation 3.12], this requirement may be written as: 

/=0 
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This polynomial in ai is equal to 0 when 



Solving a{ for pi yields the orthogonal values of /3i: 

Trk 

/?j^ = cot— , ^^ = 1,2,...,M-L (3.21) 

An interesting side observation here is that the (M - l)-element LENS parameter 
vector that consists of these orthogonal values, = [Pt.Pt^^ - ^Ph^iV^ Yi^^ds a 
LENS weight vector that is equal to the uniform-weighting vector: 



This result reinforces the choice of Qi as the 0{ robustness metric, since the white- 
noise gain of W = ^1 is = = M [Equation 2.15]. This is the maximum 
possible white-noise gain value for a weight vector that preserves ff^^rget = i with 
unit gain, which means that it achieves very robust beamforming. 

Practical Parameter Classifications 

In practice, the contours of constant are complicated and the preceding classi- 

fication procedure can be difficult to implement. In order to yield a practical system, 
this research simplifies the problem by approximating the contours as boxes, cen- 
tered at Pi = 0. Empirical observations demonstrate that boxes with the following 
proportional relationship yield reasonable matches to the desired contours: 

Height _ \Im{i}\ _ MR) ^^tsg^^'^^ 

Width \Re{P}\ R ~ R ' ^"^—^ 

where = cot ^ is the largest orthogonal /?,• value, as defined in Equation 3.21, and 
where the box half-width R lies between 0^ and oo. Based on A{R), the rectangular 
approximation becomes square as -> oo and a flat line extending from real values 
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~0t — 0M-1 to Z?/- as /? /3^. This reflects the behavior of the true contours of 
n,-, which become circular as Cli 0° and which converge to the Af - 1 points 0^, 
fc = {l,2,...,M-l},asQ,•H.90^ 

Figure 9 demonstrates the relationship between the real and approximated con- 
tours of constant for an M = 3 element array. As shown, the boxes yield 
reasonably good approximations to the actual contour curves, especially as fi, ap- 
proaches 90°. Similar results are obtained for other values of M. 

Given the box-approximation to the contours of constant fi.-, the LENS parameter 
classification from Equation 3.20 may be re-expressed as: 



l-Re{A-}| > Rt or \Im{/3i}\ > A(/2r) 
\Re{0i}\ < Rj AND \Im{Pi}\ < \{Rj) 
\Re{pi}\ > Rj OR |/m{/?.}| > \{Rj) 
AND 

\Re{(3i]\ < Rr AND |/m{A}| < A{Rr) 



A € r, 

A € J, 



(3.23) 



/?.€ru J, 



which is shown graphically in Figure 10. The values of Rr and Rj yield classification 
boundaries that are close to those obtained using ideal contours at Qr and Qj in 
Figure 8 . By definition, T and J should be non-overlapping, which is satisfied if 
Rr > Rj. In this case, the boundary-selection process becomes that of choosing Rr 
and Rj rather than that of choosing fir and ilj. These choice are explored in more 
detail below for ULAs. 



Special Case: Broadside ULA Architectures 

This section explores the process for determining the specific values of Rr and Rj 
for the special case of broadside ULA systems. 

Section 3.3.1 demonstrates that, for ULAs, there is a correspondence between the 
LENS parameters A and the beam-pattern null locations. Specifically, the parameter 
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/?,- produces a null at the location that satisfies 

= cot , (3.24) 

where r{9, (f>) is the inter-element delay experience by a source from the {6, <f>) location. 
Equation 2.7 may be used to show that ULAs that are oriented broadside to the target 
signal have inter-element delays with the structure 

Ti6,(j>) = ^sin^, (3.25) 

where 9 is the angle between the source propagation direction and the plane perpen- 
dicular to the array axis [Section 2-2]. Note that r(^,0) depends on 9 only and not 
on <t), due to the cylindrically-symmetric array behavior about the array axis. 

Equations 3.24 and 3.25 provide a convenient means for determining the target- 
canceling and jammer-canceling region boundaries that are used in the LENS classi- 
fication procedure above. Specifically, a that produces a null within an angle 9r 
of the target location is labeled target-canceling, while a /3i that produces a null at 
an angle greater than 9j from the target location is labeled jammer-canceling. Given 
that the target location is Ox = 0, these requirements may be expressed as: 

Target-Canceling: - 9t\ < 9r ^ |AI > cot sin 9r) , 
Jammer-Canceling: \9i - 9t\ > 9j |^,| < cot (gsin^j-) . 

By comparing this classification to that proposed in Equation 3.23 and by noting that 
Pi are assumed to be real for Equation 3.24, then this argument suggests that the 
desired robustness regions could be achieved using: 

Rr = cot it siner) 
Rj = cot [tsm9j). 

As long as 9r < 9j, these values will yield three distinct LENS classification regions 
as shown in Figure 10 and described in Equation 3.23. This research uses these 
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formulae for Rr and Rj, with Or = IS"" and 6j = 30^ 

In practice, both upper and lower limits are imposed upon Rr and Rj, The 

upper limit restricts any unusually large (5i values, which prevents very large values 
of the weight vector norm, EiLENS,A/(^)^J3£LENS,A/(^) prevents excessively small 

values of the system white noise gain = Ty 77^ [Equation 2.15], 

It must be noted that the precise relationship between 0 and the weight vector norm 
is complex, and so the upper bounds on Rr and Rj are not designed to restrict 
li^LENS^Af (^)^1ElenSwV/(^) be below a specific value . This research uses Rr < 30 
and i?jr < 10, which have been determined empirically to yield good performance. 

The lower limit on Rr and Rj ensures that the LENS parameters always has 
a minimal amount of 'maneuvering room' about the orthogonal parameter values 
0^ = ^^^^-^ k = {1,2, ...,M - 1}. The precise impact of these lower bounds is 
difficult to derive, and this research uses Rr > 0^ + 2 ^nd Rj > 0^ + 1, which have 
been determined empirically to yield good performance. 

Together, the upper and lower bounds on Rr and Rj for this research are: 

(3t + ^<RT< 30 

(3.27) 

Pt- + l<Rj< 10. 
3,3.3 Final Transformation: 0 ^ (3 ^ 

C_opt C—roD 

Given the LENS classification procedure above, it is possible to state the actual 
i^pt ^ob transformatiofx, based on Equation 3.11: 



0rob,i — < 



0 i3opt,i € r 

(^opt,i 0opt,i € J (3.28) 
F[/3,pt,i] 0opt,i € TUJ 



The first two components of this transformation are straightforward given the 
classes T and J. For the third component, recall that Fy3opt,i] should be a smooth 
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10 



15 
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transition between the /3opt,i € T and 0opt,i € J extremes. For the box classifications 
of Equation 3.23 [Figure 3-7], one possible F[/3opt,i] is: 



= min 



\Re3o,u\ - Rj \Im0opt.i\-A{Rj) 
Rr-Rj ' A{Rr)-HRj) 



opt, it 



(3,29) 



where f[/3optA = 1 for Popt.i on the border of the J, /[/?opt,:] = 0 for Popui on the 
border of the T, and flPopu] linearly moves from 1 to 0 as either Re[/3opt^i] progresses 
from Rj to Rr or lin[(5optA progresses from K{Rj) to !s.{Rt). Given this F[/?opt,:], 
Equation 3.2S may also be expressed as 



0robs = f[PoptAPopi,u where f[3opt^i] = < 



3.3.4 Example Cases 



0 
1 

fiPoptA 



f^opui e r 

0opt,i € 



(3.30) 



In order to gain some understanding into how the LENS robustness restriction actu- 
ally operates, this section presents some simple examples based upon a two-element 
array with an element spacing of c? = 0.07 m. This array is oriented broadside to the 
target location, which means that no steering delays are needed to make target = i- 
Assuming far-field sources, then Equations 2.6 and 2.7 may be combined to yield the 
source-to-array transfer function expression: 



H{9, 4>) = 



1 




gju/ ^ sin $ 





0.07 m n 
345m/«c ^'"^ 



where an acoustic propagation velocity of c = 345 m/sec has been used and where 6 is 
the angle between the source propagation direction and the plane that is perpendicular 
to the array axis [Section 2.2.1]. In the azimuth plane, 9 may be interpreted as the 
angle to the left of the assumed target location. 

Given this array, consider the following four source situations, where a single 
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source, with source-to-array transfer function , dominates ^^pui- The values of 
/3^p^ 1 that yield W_hEm^i{(^opui)^ = 0 ^Iso given, as determined from Equa- 
tion 3.3. Note that a perfectly-aligned target source, with Hi9^4>) ~ target? not 
considered, since, by definition, I£lENS,2(A) cannot null this location. 

1. / — 500 Hz, 10'' mis-aligned target, matched array elements: 



j27r50OHz ,°°^"^ sin 10* 

g** 345 m/ sec 



0Ui = is.oi. 



2. / = 500 Hz, aligned target, mis-matched array elements: 



0.9 



3. / = 500 Hz, 45° jammer, matched array elements: 



1 

pi2-50OHz5^f^si„45<' 



/?opM = 4.36. 



4. / = 500 Hz. 45° jammer, mis-matched array elements: 



Q^girc600Hz55°tf7^si„45<' 



/3,Vi = 4.13-il.O. 



Figure 11 plots these four 0^^^ ^ values as well as the LENS robustness regions for 
/ = 500 Hz, defined according to the Rr and Rj of Equation 3.26 with 9r = 15° and 
9 J = 30°. 

These plots demonstrate the LENS robustness restriction at work. The target- 
canceling LENS parameters, /S^pf i and /J^p^ lie in the T region, which means that the 
robustness restriction transforms them to the neutral value of /J^p^ ^ = = 0. On 
the other hand, the jammer-canceling LENS parameters, P^^t^^ and /?opt,n He within 
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The Overall LENS Algorithm 



Figure 12 summarizes the overall LENS algorithm. 



Since this 



research will study ULAs that have been oriented broadside to the target location, 
the robustness transformation boundaries derived in Section 3.3.2 will be used: 



where d is the inter-element array spacing, c is the propagation velocity, and (5^ is 
the largest magnitude orthogonal value [Equation 3.21]. 

One LENS feature is immediately obvious when comparing LENS to traditional 
CM\' processing from Section 2.3.2: LENS separates the robustness constraining pro- 
cess from the beamforming power minimization, while CMV combines the robustness 
constraints into the beamforming power minimization. It is the separation of these 
two components that allows LENS processing to exhibit direct and flexible control 
o\*er the beamforming. 




subject to 



/3^ + 2<Rr< 30, 
(3^ + l<Rj< 10, 



(3.31) 
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5 Modified LENS: Skew LENS Processing 



Skew LENS processing is a relatively simple modification of standard LENS pro- 
cessing. Specifically, it consists of the following four steps: 

1. Given the standard LENS input, which has been steered to have fl^target = i' 
form the skewed input signal Xc = where: 



1 0 

0 C 



0 



0 C 



M-1 



(5.1) 



The matrix Z is the skewing matrix and C is the skewing parameter. 

2. Generate the optimal skewed LENS parameter set §^^^^ based on X,, i.e., use 
the standard LENS definition [Equation 3.8], with Rx^Xf^ rather than Rxx- 
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3. Restrict /? . ^ for robustness to form /?. , . 

4. Form the weight vector IilLENS,A/(^ ro6) ^^^^ form the beamformer 
output Y = ]$:LENS,A/(^,,,J''iC<. 

Given this description, it is clear that Skew LENS reduces to standard LENS pro- 
cessing when C = 1 2.nd Z becomes the identity matrix 

It is very important to note that the optimal Skew LENS parameter set, is 
defined using the standard LENS definition [Equation 3.8] applied to X^, i.e., 
is designed to preserve sources with K{&'i4>) — L Given the Skew LENS definition, 
however, it is clear that the assumed target location now has a transfer function of 



10 



^,target — 2'i — 



1 

c 

A/-1 



(5.2) 



15 



which means that ^ may actually cancel the assumed target location! As described 
below, however, the introduction of this controlled and very structured error into 
target couceutrates target-canceling beamforming information into a single LENS 
parameter, which can improve the overall LENS performance. 

The following sections present the logic behind Skew LENS processing, discuss the 
selection of the skewing parameter develop a modified LENS robustness restric- 
tion for Skew LENS pi-v^;cessing, and present complete Skew LENS description. 



5.1 Skew LENS Basics 

2Q For the special case of idealized narrowband far-field sources with a mis-aligned target. 
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LENS algorithm can completely separate the target-canceling and jammer-canceling 
beamforming effects within the optimal LENS parameter set. This allows LENS to 
prevent target cancellation without limiting jammer cancellation, which results in ar- 
ray gain levels that are nearly-equal to those observed with a perfectly-aligned target. 

For all other array environments, LENS cannot complete!}' 
separate target-canceling from jammer-canceling information in the optimal LENS 
parameter set. This means that LENS cannot restrict target cancellation without 
simultaneously limiting jammer cancellation, which prevents LENS from maintaining 
near-ideal performance levels. Skew LENS processing is an attempt to increase the 
separation of target-canceling and jammer-canceling information in the optimal LENS 
parameter set for all array environments and for all sources of target-cancellation. 

Skew LENS processing achieves this goal through the skewing matrix Z. Specif- 
ically, Z concentrates the possible target-canceling LENS parameter values into a 
relatively small subset of the overall LENS parameter space. By concentrating the 
range of possible target-canceling LENS parameter values, Skew LENS processing al- 
lows these values to be covered by a single LENS parameter. This frees the remaining 
parameters for jammer cancellation and increases the separation of target-canceling 
and jammer-canceling information within 8, 

To see how Z concentrates target-canceling LENS parameters, recall the behavior 
of these parameters. Section 3,3.1 shows that P describes a basis, {V^{§)}, for the 
null-space y[KlENSMi(l)l KLEmM§)^ ^^^re 



, i-={l,2,...,Af-l}. 



A target-canceling LENS parameter, corresponds to a V^{0) that is close to the 
Apply this observation to standard and Skew LENS processing: 
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1. For standard LENS, [target = iL this closeness occurs when 



2. For Skew LENS, target = Zl] this closeness occurs when 

f-ftTJ ^ ft -J—. 

assuming C 7^ !• 

Now, assume that the actual target-to-array transfer function, target? near to, but 
not equal to, JLarget = L For standard LENS, as target ^^ries from target? '^^ is clear 
that that the target-canceling LENS parameters will cover the entire large- region. 
For Skew LENS, on the other hand, as A:,target = -^&arget varies from ^.target, it is 
clear that the target-canceling LENS parameters will be clustered about = j^. 
This demonstrates that, as desired, Skew LENS concentrates the target-canceling 
LENS parameters. 



5.2 Selecting the Skew Parameter , C 

The Skew LENS parameter, ^, determines the amount of concentration that Z in- 
troduces into the target-canceling LENS-parameter values and the overall ability of 
Skew LENS to separate target-canceling from jammer-canceling information. In gen- 
eral, as C deviates further from the non7skewed case of ^ = 1, more skewing and 
more concentration occur. Care must be taken, however, when selecting and the 
following rules should be observed: 

1. C should not yield ZHi9,(f>) = I for any non-target (g,0): 

Since that P^^^^^ is defined using the standard LENS optimization, it preserves 
sources with source-to-array transfer functions of ZHX9,(f>) = L Given the 
beamforming goal of jammer signal attenuation, it is clearly undesirable for Z 
to map a jammer source to have ZHX9, (f>) = L 
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Consider as an example broadside ULAs. As shown in Equation 2.6, these 
arrays have non-skewed source-to-array transfer functions of the form: 



H{9, cf>) = 



1 

QjuTi{e,<ii) 



gXM-l)Ti(5,0) 



10 



15 



20 



where t{0, <?>) = f sin^, d is the inter-element spacing and c is the propagation 
velocity. Given this Kid,^), a skewing factor of C = e-^''^''(''''*'> would yield 
ZH{du<f>i) = 1, and that Skew LENS would preserve sources from the {61,61) 
location! For ULAs, ZH{e, <f>)^lis avoided when [CI 7^ h 

2. C should not skew the processing excessively far from standard LENS 
processing, i.e., |C - 1| should not be too large: 

As C deviates further from the non-skewed C = 1, the concentration of the target- 
canceling LENS parameter values increases. Unfortunately, if ^ deviates too far. 
then the skewing also begins to concentrate the jammer-canceling information 
into the same set of LENS parameter values. This reduces the separation of the 
target-canceling 3.nd jammer-canceling information. 

Consider as an example the broadside ULA situation from above, which has the 
non-skewed source-to-array transfer functions given above. At low frequencies, , 
the delay term t{$, 4>) is small for both target and jammer locations, which 
means that (f>) is always 'near' to target = L A skewing parameter that 
deviates largely from C = 1 masks entirely the eP'^'^i^^^'^ term in ZHiO, <t>), which 
means that the jammer-canceling LENS parameter values for this location are 
concentrated together with the target-canceling LENS parameter values. Em- 
pirical observations conducted for this research suggest that < 0.05 avoids 
excessive concentration of jammer-canceling LENS parameter values. 



25 For this research, Skew LENS is implemented with C = 0.95. 
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5.3 Skew LENS Robustness Restriction 



As described above, Skew LENS processing applies the skewing matrix Z to the 
standard LENS beamformer input, X,, in order to concentrate the possible target- 
canceling LENS parameter values into a relatively small subset of the LENS parameter 
space. Since this transformation re-maps the target-canceling and jammer-canceling 
LENS parameter regions, it is necessary to re-define the LENS robustness regions T, 
J , and T\J J (target-canceling, jammer-canceling, and intermediate, respectively) 
that are used in the LENS robustness restriction of Section 3.3. 

Such a re-definition would be a time-consuming and rather difficult procedure, 
but, fortunately, there is a relatively simple means of by-passing this problem. It is 
possible to convert into the non-skewed LENS parameter set ^^^p^, for which 
the target-canceling and jammer-canceling LENS parameter regions are mapped back 
to their original, non-skewed, standard-LENS locations. This 3 may then be 
restricted for robustness using the standard LENS procedure of Section 3.3 to yield 
^^^^j, and this result may be 're-skewed' to form ^ 

At first glance, it may appear that this non-skewing transformation negates all 
of the beneficial effects of Skew LENS processing and that ^ may simply equal 
the original, standard LENS parameter This behavior is not, in fact, true - all 
beneficial effects of Skew LENS processing arise out of the actual power minimization 
that yields ^^^pj, and the non-skewing transformation is simply a tool that assists in 
the robustness resiriction to form /? 



5.3.1 Relationship between and the Non-Skewed /3 



i — as 



A Skew LENS parameter set ^ and its non-skewed counterpart ^ are considered 
equivalent if they describe the same beamforming behavior. Specifically, these two 
parameter sets will be considered equivalent when the beamforming accomplished by 

I£LENS.Af (^)''2Q = I£LENS.A/(^)''2i 
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is equivalent to the beamforming accomplished by 



10 



where Z is the skewing matrix of Equation 5.1 and the relation = ZX^ has been 
used. Such beamforming equivalence is attained when 

AqZ^l£LExNS.vf(^)] = ^V[i£i,ENS,M(^3)] (5-3) 

Equation 5.3 is true when 0^ and 0^^ are related using 
the following transformation pair: 



(5.4) 



Skewed Non-Skewed ^n,,, = ji^^tl^^^^ 
Non-Skewed ^ Skewed = j^^^tm<^ 

5.3.2 Skew LENS Robustness Restriction Summary 

Figure 13 presents the overall robustness restriction process for Skew LENS pro- 
cessing. It consists of the standard LE.\S robustness restriction nested between the 
3^ 3 and 3 , ^ , transformations. 

l—i,,opt : — ns.opt ■ — ns,roo '—Q,roo 



5.4 Overall Skew LENS Algorithm 

Figure 14 shows the overall Skew LENS flow diagram. This resembles the standard 
LENS flow diagram from Section 3.4, but it adds the input-skewing step as Step 1 
and modifies the LENS robustness restriction as described by Figure 13 . 
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6 LENS Implementation 



The main computational component of either standard LENS [section 3] or Skew 
LENS [ section 5] processing is the determination of the optimal LENS parameter set, 
0^^^ or ,3^ Once this optimal parameter set has been determined* it is a straight- 
forward process to restrict the parameters for robustness and to generate the final 
beamforming weight vector. This chapter develops and simulates two implementa- 
tion methods that use the temporal stream of narrowband array input observations, 
{Xj^^^} to form a running estimate 0^^^ of the optimal LENS parameter set. Note 
that, although this chapter uses the terminology of standard LENS processing, the 
identical methods may be used to obtain 0^1^^ from {Xi!^} for Skew LENS processing. 

Recall from section 3 that the optimal LENS parameter set is defined as the 
parameter set that minimizes the expected beamformer output power that arises 
when using the LENS weight vector I£lENS,m(^): 

Lpt = . KlENSMiifRxxKj^E^SMil^)^ (6.1) 

where Rxx = E[ XX ^] is the input covariance matrix and P[P] is defined to be the 
expected beamformer output power. As shown in Section 3.2.1 [Equation 3,10], the 



closed-form solution to this optimization may be written as: 

= -roots (^^) , (6.2) 

where the operation roots(2) returns the set of roots for the polynomial with coeffi- 
cients given by the elements of the vector p. 

Since Rxx is unknown and time-varying, Equation 6.2 generally cannot be used 
to generate a single that describes the entire input. Instead, must be es- 
timated actively from the temporal stream of input observations, {X}^^}. The most 
obvious LENS implementation would use {X^^^} to formulate a time-dependent input 
covariance estimate Uj)^, which, through Equation 6.2, would yield a time-dependent 
optimal parameter set ^^J^, Unfortunately, the matrix inversion and polynomial root 
calculation render this system computationally expensive in practice. This chapter 
develops two alternative, iterative solutions to the LENS optimization of Equation 6.1 
and uses these solutions to develop actual LENS implementations. 

Relaxation Approximations 

This class of adaptive solution arises out of the Gauss-Seidel relaxation-iteration 
for linear systems of equations and solves Equation 6.1 by cycling through 
the LENS parameters and minimizing P[^] with respect to each /?,- individually 
while holding the remaining /?^, m i, fixed to their most recent \^ues. 

Gradient Approximations 

This class of adaptive solution calculates the gradient of P\fi] with respect to 
the individual j3i parameters and follows the gradient to minimize P[ff\. 

6 1 Useful FormulaP- ^^^Jg] ^'^^^ 

b.l usetul i?ormulae. , ^^^^ , —rr, and ^ 

Before developing the Relaxation and Gradient LENS implementations, this section 
presents expressions for the partial derivatives and where 

A = 0R,i + j/3/,z- These partial derivatives are required for the development of both 
iterative LENS solutions. 
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Given 



where the parameters Bi and Ci are given by: 



(6.3) 



0 

0 



(6.4) 



In this definition: (1) is the transformation matrix from Equation 3.7 that trans- 
forms the LENS beamforming weight vector W i^NSm^^ alternative beam- 
forming weight vector I[leNS,a/(^) [Appendix B] and (2) I^leNS,a/,:(^) a vector 
that is entirely independent of (3i and is defined to be: 

Af-2 

[ 1 A+i ]*-•*[ 1 ^M^x ]• (6-5) 
Note that I£leNS,a/(^) [Equation 3.6] and £^LENS,A/,i related through /?i: 



t^LENS,A/{/^)'' - \ [ t/LENS,A/.(^)'' ° ] f [ ° ^LENS.W.<^)^ ] • (6-6) 



6.2 Relaxation LENS Implementation 

Relaxation LENS processing solves for Iteratively by cycling through the LENS 
parameters one-by-one and minimizing Equation 6 J with respect to each individ- 
ually, while holding the remaining parameters fixed. This process is based upon the 
idea of the Gauss-Seidel iteration to solve linear systems of equations. 
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In order to illustrate how this relaxation iteration operates, this section begins by 
exploring Equation 6.1 in the case of a known and fixed input covariance matrix i?A'A'- 
It then generalizes the result to yield the Relaxation LENS implementations that 
combine the relaxation iteration with a time-dependent input-covariance estimate 
^A-A- to generate a time-dependent optimal LENS parameter set estimate ^J^. 

6-2.1 Covariance Relaxation 

Recall the basic LENS problem from Equation 6.1: 

£. = arg min PU31 

The relaxation LENS solution to this minimization may be stated as follows. 

1. Initialize the process with an arbitrary parameter vector such as the 



orthogonal LENS parameters §^ [Equation 3.21] . 



2. Update into '^Jk + 1] by cycling through the LENS parameters A' 

exactly once and minimizing with respect to each /?,• while holding the 
remaining (3^, m f: z, fixed at their most recent values: 

0opt.i[f^ + 1] = arg min P[(3]\ . f6 7) 

A6C ^—^l^m.>n#i fixed as desired ^ ' 

Consider this update process in detail for an arbitrary Since updates have 
already been performed for /?^, m < z, the most recent values of /?^, m i are: 

a _\ ^opt,m[k +1] m < 2, 

[ Popt.mlk] m > i. 

Given these fixed f3^,m^ i, it is possible to form 5.- and d [Equation 6.4] and 
to use these to define |^ and ^ [Equation 6.3]. Since 
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then the solution to Equation 6.7 is obtained by setting |^ and to 0 and 
solving: 

|g = 2B„, + 2/3„,C,=0 =i. /3r„ = -^, 1 ^ .__B: 

This yields /?opt,t[fc + 1] - -ff to solve Equation 6.7, 

As oo, ^pJA;] converges to the true underlying 0^^^ [as demonstrated by the 

simulations in Section 6.5.1], Figure 15 depicts the overall LENS relaxation iteration 
process. 



6.2.2 Relaxation LENS Processing 

10 Now that the basic concept of the LENS relaxation iteration in the case of a known 

and fixed input covariance has been derived, it is necessary to modify the iteration 
to suit the more realistic situation in which the input covariance is unknown and 
time- varying. This is achieved by using the temporal stream of array inputs, {X^^^}, 
to form a time-dependent estimate R^j^j^ of Rxx^ 

15 The most direct Relaxation LENS implementation would use each r!xx ^^rm 

a corresponding P[y3]<^) and would apply a full relaxation solution until ^^j^ has been 
generated to the desired precision. Assuming that the relaxation iteration converges 
to the true solution to Equation 6.1 (this is demonstrated below), this system is 
equivalent to one that uses Equation 6.2 directly to generate from -RxA'- Unfor- 

20 tunately, this implementation has the disadvantage that it requires several iterations 

to obtain each which may be computationally demanding. 

Simplified adaptive relaxation implementations result by realizing that r!xx 
correlated with RIxx^*^ many of the same input observa- 

tions go into the averages that form e!^x ^xx^^ ^iiis reason, and ^^J'^ 

25 are generally 'close' to one another. The number of relaxation iterations required to 

obtain can be substantially reduced by using as the initialization for the 

iteration to solve for In fact, if ^'^^ is sufficiently slowly-changing, then it is 
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possible to update only a fraction of the LENS parameters at each input step while 
maintaining reasonable solution accuracy. 

Figure 16 shows two possible simplified implementations, using the notation 
^MPopt.,] to denote a relaxation iteration update of the parameter ^^J/,,,-. 

^- '^^^ Complete Update system updates each of the M - 1 LENS parameters in 
^ . to yield ^ . 

2. The Partial Update system updates exactly one of the LENS parameters in 
to yield ^'J^. Note that for this system, care must be taken so that the 
relaxation update still cycles through all M-1 LENS parameters. This cycling 
is achieved by updating parameter z' = [1 + (/ mod M - 1)] at input step /. 

The primary trade-off between these two implementations is one of implementation 
complexity versus solution accuracy and adaptation speed. By updating only one pa- 
rameter per input observation, the partial update implementation reduces the compu- 
tations required by a factor of relative to the complete update implementation. 
This reduced complexity comes at the cost of updating each LENS parameter only 
once every M - 1 input observations, which increases the system adaptation time to 
a changing environment and reduces the overall solution accuracy, .^s shown by the 
simulations of Section 6.5.2, however, the co\'ariance estimate R$j. generally changes 
slowly enough that this latter effect is minor and that Complete and Partial Update 
Relaxation LENS implementations yield nearly identical results. 



6.3 Gradient LENS Implementation 

Gradient LENS processing takes a different approach to solving Equation 6.1. Specif- 
ically, it calculates the gradient of P[0} with respect to the LENS parameters and 
steps along this gradient until P[0] has been minimized. In a manner similar to that 
of the previous section, this section develops gradient LENS processing by consider- 
ing the problem in the case of a known co\'ariance matrix, and then generalizing the 
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result to yield the Gradient LENS implementation that generates the time-dependent 
0ipt for an unknown covariance matrix. 



6.3.1 Known Covariance Gradient Searcii 

Recall the basic LENS parameter problem from Equation 6,1, 



0 = arg min PW] 



The basic gradient-based LENS solution to this minimization may be stated as follows. 

1. Initialize ^^JO] with an arbitrary parameter vector, such as the orthogonal 
LENS parameter set, from Equation 3.2L 

2. Update p^Jk] into P^Jk + l] by calculating the gradient of P[p] at p^Jk] and 
taking a small step along this gradient toward the function minimum: 



[k] ' 

—opt*- J 



(6.S) 



where Equation 6.3 allows VP[,/3] to be written as: 



VP[0\ = 


dPW] 

80R,2 


+ j 


' dPl0] ■ 

aPi/3] 

301.2 


= 2 






om 




dP[0] 












. d0I,M-i J 





Using the definitions of B,- and d [Equation 6.4], the vectors C/;LENS.iw(^) and 
LLle:^SmA§) [Equations 3.6 and 6.5], and the relationship between £lenS, A/ {^) 
and Uiens,mA§) [Equation 6.6], B' + f3iC, may be re-written as: 



£lens.w,;(^)^ 
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where the second equality arises from the fact that {B* + (3iCi)^ = S* + Ad- 
This expression for 5* + pid yields: 



0 (LlEmMAlf 
0 £lLENS.Vf,2{^)^ 



^XT;fT^LENS,Af(£)*- 



(6.9) 



The step-size parameter {i in the gradient update [Equation 6.8] governs the over- 
all rate and quality of the iteration convergence. A larger step-size results in faster 
initial convergence towards the optimum, but exhibits oscillation about the optimal 
value (also known as misadjustment). A smaller step-size results in slower conver- 
gence but less misadjustment. It should also be noted that, if the stejp-size is too large, 
then this iteration may diverge. This research presents no proofs nor performs any 
investigations into the range of step-size that guarantees convergence of this LENS 
gradient-search. Simulations suggest, however, that the general gradient-based con- 
vergence rule 

Iteration a 
Convergence 

yields acceptable results provided that the parameter magnitudes, |/?opt,j[^jl » i = 
{1, 2, • ■ • , M — 1}, are not too large. This research imposes the following modification 
upon the step-size in order to reduce the step size when the |j5opt,t[fc]P large: 

= A'sd^^ - 0<a<I, 
with A- = max[l, ^.= - /_7-W'' ]. 

In this case, the parameter A' is restricted to be greater than 1 to ensure that it only 
shrinks and never enlarges fi[k]. 
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6.3.2 Gradient LENS Processing 



The most direct gradient-search modification to accommodate an unknown and chang- 
ing input covariance would use the input stream, {2L^^^} to form R^^^ 

to apply a full gradient search to each ^Ja' V'^^^^ Such a 
system would be inefHcient, however, since gradient search iterations generally have 
long convergence times^ 

A simpler implementation may be realized by synchronizing the gradient-search 
iteration to the input observation iteration. Specifically, use r!xx generate the 
current gradient VP[^(') [Equation 6.9] and substitute this into Equation 6.S to 
update f\ into f^^^^: 

€"=i!p.-''"'^^ta"'U... (6.11) 

For this implementation, the time-dependent step-size parameter defined in Equa- 
tion 6.10 may be approximated by: 

'^'' = ^ 0<a<l, 
with A' = maxfl, ^'='J_°^'-'^ ], 

where (T^ Ms a running estimate of the total expected array input po%ver, E[2L"}Q.K 

In general, each update step m Gradient LENS processing (or in any gradient-based 
adaptive algorithm) is smaU. Therefore, the adaptive Gradient LENS processing does 
not require a highly accurate estimate of the mput covariance to yield acceptable 
results. The Gradient LENS processing estunates igas the instantaneous mput 
covaraince: 
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(6.12) 



This reflects the covariance matrix estimate used by the traditional Frost Process- 
ing and Scaled Projection Algorithm gradient-based beamforming implementations 
[Section 2.3.3]. 

6.4 Computational Complexity 

Figure 17 plots the number of real multiplications required per input snap shot, . 
^ , by Complete Update Relaxation, Partial Update Relaxation, and Gradient LENS S 
Processing. 

Notably, the present mvention separates the robustness constraining process 
from the beamforming power mmimization, in contrast to prior art techniques which 
combine the robustness constraint into the beamforming power mmimization. One of 
ordmary skill will recognize that the signal processmg technique of the present 
invention has a wide variety of system applications. 

Although the present mvention has been shown and described with respect to 
several preferred embodunents thereof, various changes, omissions and additions to die 
form and detail tijereof, may be made therein, witiiout departing from the spirit and 
scope of the invention. 

What is claimed is: 
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CLAIMS 



11. A method of selecting the tap weights W^for an adaptive multi-tap frequency domain 

2 digital filter that processes an input signal vector X from a plurality of spatially separated 

3 transducers that detect energy from a plurality of sources including a target energy source and at 

4 least one non-target energy source, wherein the filter receives and processes the input signal vector 

5 X to attenuate noise from non-target sources and provides an output signal vector Y, the method 

6 comprising the steps of: 

7 parameterizing each of the tap weights W^^ such that each of the tap weights is 

8 characterized by a vector of parameters ; 

9 solving for each parameter of the vector fi^^^ by minimizing the expected power of the array 

10 output signal Y; 

1 1 applying a robustness-control transformation to the vector to provide a robust vector 
■^^ ^rob' wherem the robustness-control transformation identifies and reduces target canceling 

13 components of the vector £,^^that arise firom incomplete target location knowledge while 

1 4 preserving non-target canceling components; and 

1 5 forming the weight vector indicative of tiie filter tap weights as a function of die vector 

16 



1 2. A signal processing apparatus that receives an input signal vector X from a plurality of 

2 spatially separated tiransducers that detect energy firom a plurality of sources includmg a target 

3 energy source and at least one non-target energy source, wherein the apparatus processes the input 

4 signal vector X with a digital filter comprising a plurality of tap weights to attenuate signal 

5 noise from non-target sources and provide a resultant output signal vector Y, said apparatus 

6 comprising: 

7 means for parameterizing each of said tap weights PF^ such that each of said tap weights 
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8 is characterized by a vector of parameters J3 : 

—opt 

9 means for solving for each parameter of the vector fi^^ hy seeking a minimum for the 

1 0 expected power of the ou^ut signal Y; 

11 means for applying a robustaess-control transformation to the vector j3 to provide a 

12 robust vector ^^^^ , wherein the robustness-control transformation identifies and reduces target 

1 3 canceling components of the vector that arise from incomplete target location knowledge while 

1 4 preserving non-target canceling components; and 

1 5 means for forming a weight vector indicative of the tap weights as a function of the vector 



16 J3 , 

—rob 
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ABSTRACT 



An adaptive multi-tap frequency domain digital filter processes an input signal vector X from a 
plurality of spatially separated transducers that detect energy from a plurality of sources including 
5 a target energy source and at least one non-target energy source. The filter receives and processes 
the input signal vector X to attenuate noise from non-target sources and provides an output signal 
vector Y. Tap weights W^^ for filter are selected by first parameterizing each of the tap weights 
such diat each of the tap weights is characterized by a vector of parameters J3 , and the 

— opt 

solving fOT each parameter of the vector by minimizing the expected power of the array output 
1 0 signal Y. A robustness-control transformation is then applied to the vector J3 ^ to provide a robust 
vector ^^^^ , wherein the robustness-control transformation identifies and reduces target canceling 
components of the vector while preserving non-target canceling components. Finally, the 
weight vector indicative of the filter tap weights is formed as a function of the vector 3 . 

drob 

Notably, the present invention separates the robustness constraining process from the beamforming 
1 5 power minimization, in contrast to prior art techniques which combine the robustness constraint 
into the beamforming power minimization. The present invention uses a direct and flexible 
robustness control mechanism to yield a beamformer that provides good performance and is robust 
to a wide variety of adverse conditions. 
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Contours for = 30° (- - -) 60° ( — ), and 90" (x) for 
M = 2,3,4, and 5 element arrays. 
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: Contours ( ) and approximated contours ( ) for M=3 

and = 30°, 40°, 50°, and 60°, with /?/• also indicated (x). 
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Diagram showing the three modified LENS parameter clas- 
sification regions for an M=:3 element array as according to Equa- 
tion 3.23 when Rr and Rj are chosen to reflect the Or = 30*" and 
fly = eO"" contours, respectively. 
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Four example values of /?opt,i for a two-element array along 
with the corresponding LENS robustness regions. 



X 

Steered so that \ 

^target = 1 / 

1 

Step 1: Solve for [Eq. 3.8] 
Step .2: Make robust [Eqs. 3.30, 3.29] 

Prob,i = f[/3opt,i]l3opt,i, where 
f 0 Popt,i€T 

JiPoptA - mm [ r^ZrJ-^^ \{Rt)-\{Rs) J 



Step 3: Form weight vector [Eq. 3.1] 

\ 

W for beamforming 



1 



pt 



Step 1: Transform to [Eq. 5.4] 

o _ • (/gC,opt.»+j)+<(/?f .opt.. +3) 

Pns,opt,i - 3 (/J<,„p,,,+j)-<{/3f.op...+i) . 



Step 2: Standard LENS Robustness 
Restriction [Eqs. 3.30, 3.29] 

= f[/3ns,opt,i]i^opt,i, where 



r 0 

fy3ns,opt,i] = < 1 



er 

0opt,i € J 



/LPns,opt.«J — [ Rr-Rj 



[ f[^ns,optA 0ns,opt,i€TUj 

\R<:i0ns.ovt.i}\-Rj |rni{i3ns.opt..}|-A(fij) 1 
' \{Rt)-HRj] J 



Step 3: Transform to ^^^^ [Eq. 5.4] 



a _ • <l0ns.robi.i+3)+{0ns.Tob.i+i) 



^,rob 



^c. 13 



X 

(Steered so that \ 
^axget = i / 
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Step 1: Skew Input [Eq. .5.1] 
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Step 2: Solve for ^^^^^ [Eq. 3.8] 

I 



Step 3: Make §^ robust 
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Step 4: Form weight vector [Eq. 3.1] 

m = ]£LENS.Af(^.,^) 

T 



W for beamforming 



for i 


= 1 : M - 1, 


1. 


Fix /?m, m 7^ f at 








\ $opt,m[k\ m> i ' 


2. 


Form £LLENS,Af,i(^) fro"^ /?m. 


3. 


Use £^LENS,A/,:(^)i ^^A', and Tm 




to form Bi and d. 


4. 


KtAk + 1] = -f . 


end 
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Form 



for z = 1 : M — 1, 
end 



£-opt 
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... ,. Input 
Initialize Stream 
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Form 
6(0 



1. i = l + ^mod(M-l) 



2. = RJx[4S,d 



(a) Complete Update 



(b) Partial Update 
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